/*
 * This file is subject to the terms of the GFX License. If a copy of
 * the license was not distributed with this file, you can obtain one at:
 *
 *              http://ugfx.org/license.html
 */

#include "gfx.h"

#if GFX_USE_GMISC

#if GMISC_NEED_FASTTRIG
	const double sintabledouble[] = {
		0.000000, 0.017452, 0.034899, 0.052336, 0.069756, 0.087156, 0.104528, 0.121869,
		0.139173, 0.156434, 0.173648, 0.190809, 0.207912, 0.224951, 0.241922, 0.258819,
		0.275637, 0.292372, 0.309017, 0.325568, 0.342020, 0.358368, 0.374607, 0.390731,
		0.406737, 0.422618, 0.438371, 0.453990, 0.469472, 0.484810, 0.500000, 0.515038,
		0.529919, 0.544639, 0.559193, 0.573576, 0.587785, 0.601815, 0.615661, 0.629320,
		0.642788, 0.656059, 0.669131, 0.681998, 0.694658, 0.707107, 0.719340, 0.731354,
		0.743145, 0.754710, 0.766044, 0.777146, 0.788011, 0.798636, 0.809017, 0.819152,
		0.829038, 0.838671, 0.848048, 0.857167, 0.866025, 0.874620, 0.882948, 0.891007,
		0.898794, 0.906308, 0.913545, 0.920505, 0.927184, 0.933580, 0.939693, 0.945519,
		0.951057, 0.956305, 0.961262, 0.965926, 0.970296, 0.974370, 0.978148, 0.981627,
		0.984808, 0.987688, 0.990268, 0.992546, 0.994522, 0.996195, 0.997564, 0.998630,
		0.999391, 0.999848, 1.000000, 0.999848, 0.999391, 0.998630, 0.997564, 0.996195,
		0.994522, 0.992546, 0.990268, 0.987688, 0.984808, 0.981627, 0.978148, 0.974370,
		0.970296, 0.965926, 0.961262, 0.956305, 0.951057, 0.945519, 0.939693, 0.933580,
		0.927184, 0.920505, 0.913545, 0.906308, 0.898794, 0.891007, 0.882948, 0.874620,
		0.866025, 0.857167, 0.848048, 0.838671, 0.829038, 0.819152, 0.809017, 0.798636,
		0.788011, 0.777146, 0.766044, 0.754710, 0.743145, 0.731354, 0.719340, 0.707107,
		0.694658, 0.681998, 0.669131, 0.656059, 0.642788, 0.629320, 0.615661, 0.601815,
		0.587785, 0.573576, 0.559193, 0.544639, 0.529919, 0.515038, 0.500000, 0.484810,
		0.469472, 0.453990, 0.438371, 0.422618, 0.406737, 0.390731, 0.374607, 0.358368,
		0.342020, 0.325568, 0.309017, 0.292372, 0.275637, 0.258819, 0.241922, 0.224951,
		0.207912, 0.190809, 0.173648, 0.156434, 0.139173, 0.121869, 0.104528, 0.087156,
		0.069756, 0.052336, 0.034899, 0.017452, 0.000000, -0.017452, -0.034899, -0.052336,
		-0.069756, -0.087156, -0.104528, -0.121869, -0.139173, -0.156434, -0.173648, -0.190809,
		-0.207912, -0.224951, -0.241922, -0.258819, -0.275637, -0.292372, -0.309017, -0.325568,
		-0.342020, -0.358368, -0.374607, -0.390731, -0.406737, -0.422618, -0.438371, -0.453990,
		-0.469472, -0.484810, -0.500000, -0.515038, -0.529919, -0.544639, -0.559193, -0.573576,
		-0.587785, -0.601815, -0.615661, -0.629320, -0.642788, -0.656059, -0.669131, -0.681998,
		-0.694658, -0.707107, -0.719340, -0.731354, -0.743145, -0.754710, -0.766044, -0.777146,
		-0.788011, -0.798636, -0.809017, -0.819152, -0.829038, -0.838671, -0.848048, -0.857167,
		-0.866025, -0.874620, -0.882948, -0.891007, -0.898794, -0.906308, -0.913545, -0.920505,
		-0.927184, -0.933580, -0.939693, -0.945519, -0.951057, -0.956305, -0.961262, -0.965926,
		-0.970296, -0.974370, -0.978148, -0.981627, -0.984808, -0.987688, -0.990268, -0.992546,
		-0.994522, -0.996195, -0.997564, -0.998630, -0.999391, -0.999848, -1.000000, -0.999848,
		-0.999391, -0.998630, -0.997564, -0.996195, -0.994522, -0.992546, -0.990268, -0.987688,
		-0.984808, -0.981627, -0.978148, -0.974370, -0.970296, -0.965926, -0.961262, -0.956305,
		-0.951057, -0.945519, -0.939693, -0.933580, -0.927184, -0.920505, -0.913545, -0.906308,
		-0.898794, -0.891007, -0.882948, -0.874620, -0.866025, -0.857167, -0.848048, -0.838671,
		-0.829038, -0.819152, -0.809017, -0.798636, -0.788011, -0.777146, -0.766044, -0.754710,
		-0.743145, -0.731354, -0.719340, -0.707107, -0.694658, -0.681998, -0.669131, -0.656059,
		-0.642788, -0.629320, -0.615661, -0.601815, -0.587785, -0.573576, -0.559193, -0.544639,
		-0.529919, -0.515038, -0.500000, -0.484810, -0.469472, -0.453990, -0.438371, -0.422618,
		-0.406737, -0.390731, -0.374607, -0.358368, -0.342020, -0.325568, -0.309017, -0.292372,
		-0.275637, -0.258819, -0.241922, -0.224951, -0.207912, -0.190809, -0.173648, -0.156434,
		-0.139173, -0.121869, -0.104528, -0.087156, -0.069756, -0.052336, -0.034899, -0.017452
		};

	double fsin(int degrees) {
		if (degrees < 0)
			degrees -= (degrees/360-1)*360;
		else if (degrees >= 360)
			degrees %= 360;
		return sintabledouble[degrees];
	}

	double fcos(int degrees) {
		return fsin(degrees+90);
	}

#endif

#if GMISC_NEED_FIXEDTRIG
	const fixed sintablefixed[] = {
		0, 1143, 2287, 3429, 4571, 5711, 6850, 7986,
		9120, 10252, 11380, 12504, 13625, 14742, 15854, 16961,
		18064, 19160, 20251, 21336, 22414, 23486, 24550, 25606,
		26655, 27696, 28729, 29752, 30767, 31772, 32767, 33753,
		34728, 35693, 36647, 37589, 38521, 39440, 40347, 41243,
		42125, 42995, 43852, 44695, 45525, 46340, 47142, 47929,
		48702, 49460, 50203, 50931, 51643, 52339, 53019, 53683,
		54331, 54963, 55577, 56175, 56755, 57319, 57864, 58393,
		58903, 59395, 59870, 60326, 60763, 61183, 61583, 61965,
		62328, 62672, 62997, 63302, 63589, 63856, 64103, 64331,
		64540, 64729, 64898, 65047, 65176, 65286, 65376, 65446,
		65496, 65526, 65536, 65526, 65496, 65446, 65376, 65286,
		65176, 65047, 64898, 64729, 64540, 64331, 64103, 63856,
		63589, 63302, 62997, 62672, 62328, 61965, 61583, 61183,
		60763, 60326, 59870, 59395, 58903, 58393, 57864, 57319,
		56755, 56175, 55577, 54963, 54331, 53683, 53019, 52339,
		51643, 50931, 50203, 49460, 48702, 47929, 47142, 46340,
		45525, 44695, 43852, 42995, 42125, 41243, 40347, 39440,
		38521, 37589, 36647, 35693, 34728, 33753, 32767, 31772,
		30767, 29752, 28729, 27696, 26655, 25606, 24550, 23486,
		22414, 21336, 20251, 19160, 18064, 16961, 15854, 14742,
		13625, 12504, 11380, 10252, 9120, 7986, 6850, 5711,
		4571, 3429, 2287, 1143, 0, -1143, -2287, -3429,
		-4571, -5711, -6850, -7986, -9120, -10252, -11380, -12504,
		-13625, -14742, -15854, -16961, -18064, -19160, -20251, -21336,
		-22414, -23486, -24550, -25606, -26655, -27696, -28729, -29752,
		-30767, -31772, -32768, -33753, -34728, -35693, -36647, -37589,
		-38521, -39440, -40347, -41243, -42125, -42995, -43852, -44695,
		-45525, -46340, -47142, -47929, -48702, -49460, -50203, -50931,
		-51643, -52339, -53019, -53683, -54331, -54963, -55577, -56175,
		-56755, -57319, -57864, -58393, -58903, -59395, -59870, -60326,
		-60763, -61183, -61583, -61965, -62328, -62672, -62997, -63302,
		-63589, -63856, -64103, -64331, -64540, -64729, -64898, -65047,
		-65176, -65286, -65376, -65446, -65496, -65526, -65536, -65526,
		-65496, -65446, -65376, -65286, -65176, -65047, -64898, -64729,
		-64540, -64331, -64103, -63856, -63589, -63302, -62997, -62672,
		-62328, -61965, -61583, -61183, -60763, -60326, -59870, -59395,
		-58903, -58393, -57864, -57319, -56755, -56175, -55577, -54963,
		-54331, -53683, -53019, -52339, -51643, -50931, -50203, -49460,
		-48702, -47929, -47142, -46340, -45525, -44695, -43852, -42995,
		-42125, -41243, -40347, -39440, -38521, -37589, -36647, -35693,
		-34728, -33753, -32768, -31772, -30767, -29752, -28729, -27696,
		-26655, -25606, -24550, -23486, -22414, -21336, -20251, -19160,
		-18064, -16961, -15854, -14742, -13625, -12504, -11380, -10252,
		-9120, -7986, -6850, -5711, -4571, -3429, -2287, -1143
		};

	fixed ffsin(int degrees) {
		if (degrees < 0)
			degrees -= (degrees/360-1)*360;
		else if (degrees >= 360)
			degrees %= 360;
		return sintablefixed[degrees];
	}

	fixed ffcos(int degrees) {
		return ffsin(degrees+90);
	}

#endif

#if GMISC_NEED_INVSQRT
	// Algorithm based on Quake code.
	#if GMISC_INVSQRT_REAL_SLOW
		#include <math.h>
		float invsqrt(float n) {
			return 1.0/sqrt(n);
		}
	#else
		float invsqrt(float n) {
			int32_t		i;
			float		x2;

			x2 = n * 0.5F;

			// Convert into an int32 (no binary format conversion)
			#if GMISC_INVSQRT_MIXED_ENDIAN
				((char *)&i)[0] = ((char *)&n)[3];
				((char *)&i)[1] = ((char *)&n)[2];
				((char *)&i)[2] = ((char *)&n)[1];
				((char *)&i)[3] = ((char *)&n)[0];
			#else
				i  = *(int32_t *)&n;
			#endif

			// evil floating point bit level hacking
			i  = 0x5F375A86 - (i >> 1);					// The quake code used 0x5F3759DF but this is better.

			// Convert back to a float (no binary format conversion)
			#if GMISC_INVSQRT_MIXED_ENDIAN
				((char *)&n)[0] = ((char *)&i)[3];
				((char *)&n)[1] = ((char *)&i)[2];
				((char *)&n)[2] = ((char *)&i)[1];
				((char *)&n)[3] = ((char *)&i)[0];
			#else
				n  = *(float *)&i;
			#endif

			n  = n * (1.5F - (x2 * n * n));				// 1st iteration
			//n  = n * (1.5F - (x2 * n * n));			// 2nd iteration for extra precision, this can be removed
			return n;
		}
	#endif
#endif

#endif /* GFX_USE_GMISC */
